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Gross primary production (GPP) is a key terrestrial ecophysiological process that links 
atmospheric composition and vegetation processes. Study of GPP is important to global 
carbon cycles and global warming. One of the most important of these processes, plant 
photosynthesis, requires solar radiation in the 0.4-0.7 pm range (also known as 
photosynthetically active radiation or PAR), water, carbon dioxide (CO 2 ), and nutrients. 

A vegetation canopy is composed primarily of photosynthetically active vegetation (PAV) 
and non-photosvnthetic vegetation (NPV; e.g., senescent foliage, branches and stems). A 
green leaf is composed of chlorophyll and various proportions of nonphotosynthetic 
components (e.g., other pigments in the leaf, primary/secondary/ tertiary veins, and cell 
walls). The fraction of PAR absorbed by whole vegetation canopy (FAPAR can0 py) has 
been widely used in satellite-based Production Efficiency Models to estimate GPP (as a 
product of FAPARcanopy x PAR x LUEcanopy, where LUE canopy is light use efficiency at 
canopy level). Elowever, only the PAR absorbed by chlorophyll (a product of FAPAR C hi x 
PAR) is used for photosynthesis. Therefore, remote sensing driven biogeochemical 
models that use FAPARchi in estimating GPP (as a product of FAPAR C hi x PAR x LUE C hi) 
are more likely to be consistent with plant photosynthesis processes. 

Our paper has been designed to test which group ([FAPARcanopy, LUE can opy] vs. 
[FAPARchi, LUEchi]) is more consistent with plant photosynthesis processes. 




Using a coupled canopy-leaf radiative transfer model, we have estimated 
FAPARchi and FAPARcanopy for the Southern Old Aspen forest (SOA) in Canada for 
2001-2005 with MODIS images. The tower fluxes over the SOA site provide real time 
photosynthesis of the forest. The scientists of the SOA site offer the measurements of 
photosynthesis and their flux tower based LUE (LUE t0W er)- Our results showed that 
LUE C hi matched well with LUE towe r both at magnitude and at phase while LUE can opy did 
not. Using FAPAR can opy to estimate absorbed PAR for photosynthesis will significantly 
overestimate. One can’t get good estimate of GPP phenology if using the group of 
FAPARcanopy and LUE can opy to estimate GPP. Using FAPAR C hi and LUE C ht to estimate 
GPP will get more consistent results and will help the study of global carbon cycles and 
global w'arming. 
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2 Abstract 

3 We used daily MODerate resolution Imaging Spectroradiometer (MODIS) imagery obtained 

4 over a five-year period to analyze the seasonal and inter- annual variability of the fraction of 

5 absorbed photosynthetically active radiation (FAPAR) and photosynthetic light use efficiency 

6 (LUE) for the Southern Old Aspen (SOA) flux tower site located near the southern limit of the 

7 boreal forest in Saskatchewan, Canada. To obtain the spectral characteristics of a standardized 

8 land area to compare with tower measurements, we scaled up the nominal 500 m MODIS 

9 products to a 2.5 km x 2.5 km area (5x5 MODIS 500 m grid cells). We then used the scaled-up 

10 MODIS products in a coupled canopy-leaf radiative transfer model, PROSAIL-2, to estimate the 

1 1 fraction of absorbed photosynthetically active radiation (APAR) by the part of the canopy 

1 2 dominated by chlorophyll (FAPAR ch i) versus that by the whole canopy (FAPAR canopy ). Using 

13 the additional information provided by flux tower-based measurements of gross ecosystem 

14 production (GEP) and incident PAR, we determined 90-minute averages for APAR and LUE 

15 (slope of GEP: APAR) for both the physiologically active foliage (APAR ch i, LUE C i,i) and for the 

16 entire canopy (APAR can0 py, LUE can o P y)- 

17 The flux tower measurements of GEP were strongly related to the MODIS-derived 

1 8 estimates of APAR C hi (r 2 = 0.78) but only weakly related to APAR canopv (r 2 = 0.33). Gross LUE 

19 between 2001 and 2005 for LUE ch i was 0.0241 pmol C jxmol 1 PPFD whereas LUEcanopy was 

20 36% lower. Time series of the 5-year normalized difference vegetation index (NDVI) were used 

21 to estimate the average length of the core growing season as days of year 152 - 259. Inter-annual 

22 variability in the core growing season LUE C hi (jimol C jimol 1 PPFD) ranged from 0.0225 in 

23 2003 to 0.0310 in 2004. The five-year time series of LUE C hi corresponded well with both the 

24 seasonal phase and amplitude of LUE from the tower measurements but this was not the case for 
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1 LUE canop y. We conclude that LUE c hi derived from MODIS observations could provide a more 

2 physiologically realistic parameter than the more commonly used LUE cas ,opy as an input to large- 

3 scale photosynthesis models. 

4 

5 Key W'ords: MODIS; aspen; chlorophyll; ecosystem flux; gross primary production; gross 

6 ecosystem production; light use efficiency (LUE); LUE C hi; FAPAR ch i; NDVI; LUE toW er 
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1 . 


Introduction 


2 1.1 Background - Using Light Use Efficiency to Estimate Ecosystem Photosynthesis 

3 Realistic models of plant canopy photosynthesis are necessary for obtaining accurate estimates of 

4 the carbon cycle for use in land surface models (LSMs) and atmospheric general circulation 

5 models (GCMs) (Sellers et al. 1996a, 1996b). In vegetative canopies, photosynthetically active 

6 radiation (PAR) is absorbed from sunlight by photosynthetic pigments, primarily chlorophyll a 

7 and its accessory pigments (chlorophyll b, carotenoids). When ecosystem photosynthesis is 

8 calculated with a process model, it is referred to as Gross Primary Production (GPP). When it is 

9 calculated from flux tower data, it is referred to as Gross Ecosystem Production, designated here 

10 as GEPtower- 

1 1 Plant production efficiency models (PEMs) have been developed to estimate GPP at 

12 canopy, landscape, regional and global scales, utilizing optical remote sensing to provide the 

13 fraction of absorbed PAR (FAPAR). Examples include GLO-PEM (Prince et al. 1995, 2000; 

14 Prince and Goward 1995, 1996), TURC (Ruimy et al. 1994, 1996a, 1996b), 3-PG (Landsberg 

15 and Waring 1997; Law et al. 2000) and PSN (Running et al. 1994, 1999a, 1999b, 2000, 2004). 

16 This latter model is a satellite-based global photosynthesis product derived from the MODerate 

1 7 resolution Imaging Spectroradiometer (MODIS) on the Terra and Aqua platforms. 

18 All of these models estimate GPP as the product of three terms: the light use efficiency of 

19 the canopy (LUE ca nopy). which is a measure of the PAR conversion efficiency into 

20 photosynthetically fixed C0 2 ; the FAPAR of the canopy (FAPAR eanopy ), which is estimated 

21 using radiative transfer models and remote sensing data or using empirical relationship between 

22 FAPAReanopy and the normalized difference vegetation index (NDVI, Tucker 1979); and the 

23 incident PAR where: 
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GPP = LUE canopy * FAPAR cannpy * PAR. (1) 

Consequently, accurate estimates of FAPAR and LUE for ecosystems are essential for obtaining 
accurate GPP. 

The LUE concept was initially developed for agricultural crops at harvest to determine 
the conversion efficiency of available light into biomass (g C dry mass) over a full growing 
season and is typically expressed in units such as g C MJ 1 PAR (Monteith 1972, 1977). This 
seasonal crop-level LUE represents a direct measure of the average conversion efficiency of all 
above ground plant material (Gower et al. 1999), which is dominated by foliage for agricultural 
crops. Eddy covariance flux towers have the capability to provide near-continuous 
measurements of GEP - denoted as GEP to wer, and absorbed PAR - denoted as APAR t0W er (see 
section 3.2.3 for more details, also see Krishnan et al., 2006), for an entire ecosystem for time 
periods as short as 30 minutes. Consequently, these instrumented flux towers also provide near- 
continuous measurements of LUE, denoted as LUE towe r, over these same time periods as: 


LUE,, 


GEP,, 


GEP,„ 


( 2 ) 


FAPAR wmr * PAR tower APAR lower 

The LUEt 0W er is typically expressed as |Ltmol C0 2 nmol' 1 PAR or jimol C nmol ' PPFD, where 
PPFD is the photosynthetic photon flux density. For these tower-based calculations, incident 
PAR is measured directly by radiometers attached to the flux tower and the FAPAR estimate is 
based both on detailed canopy structural measurements and on the radiometer measurements 
(Barr et al. 2007). An underlying assumption supporting the LUE retrieval through the MODIS 
modeling approach is that the LUE can0 py used in the models is a good approximation of LUEtower, 
at least when the measurement footprint of the tower is roughly equivalent to the area of the 
remote sensing pixel. Apparent ecosystem LUE obtained from flux tower measurements (i.e., 
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1 LUEtower) directly describes the integrated response of the whole ecosystem to the prevailing 

2 environmental conditions, as do remotely acquired spectral snapshots although these latter are 

3 limited to specific acquisition times and viewing configurations. 

4 On a canopy or ecosystem scale, GEP and APAR are typically linearly related (e.g., 

5 Waring et ah, 1995), so that LUE can be determined from the slope of this relationship. This 

6 apparent linearity results from multiple scattering within the canopy, w'hich involves 3-D 

7 contributions from foliage of multiple species and illumination conditions, as well as non- 

8 photosynthetic material (e.g., limbs, trunks, cones, litter). On the other hand, comparable light 

9 response curves for individual leaves of selected species yield non-linear responses for which the 

10 initial slope of the linear portion of the curve describes the quantum efficiency (Mohr et al., 

1 1 1995). The quantum efficiency of individual leaves can also serve as an input to carbon cycle 

12 models but a means of scaling it to the canopy level is still required. 

13 A common modeling approach is to set a maximum LUE for optimal environmental 

14 conditions (i.e., unstressed vegetation) and to simulate ecosystem responses when unfavorable 

15 environmental conditions occur (e.g., limitations of temperature, humidity, soil moisture, etc.) 

16 through down-regulation of the maximum LUE to achieve an apparent LUE (Medlyn, 1998). 

17 The MODIS GPP product, an output of the PSN model, has been compared with 

18 measurements made at flux towers by several research groups. For instance. Turner and 

19 colleagues (Turner et al, 2003, 2004, 2006) found that the annual MODIS GPP totals calculated 

20 using MODIS standard photosynthesis products for a deciduous forest in Massachusetts, USA, 

2 1 matched well with the annual GEP totals from the flux tower. However, the seasonal time course 

22 of MODIS GPP dynamics differed significantly from the GEP measured by the flux tower 
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1 (GEP tower ) suggesting that a more physiologically realistic method of estimating GPP could be 

2 useful. 

3 

4 1.2 Chlorophyll-based LUE ( LUE C ht ) 

5 Even though maximum leaf LUE can be strongly influenced by leaf chlorophyll 

6 concentration (e.g., Waring et al. 1995), it is less clear how canopy chlorophyll concentration 

7 might influence apparent LUE at the ecosystem scale. Laboratory studies (Yoder and Waring 

8 1994) have shown that variation in canopy total chlorophyll content of miniature Douglas fir 

9 canopies was significantly correlated with their photosynthesis, although the correlation was 

10 higher for canopies exposed to full sun. Several other studies have shown a relationship between 

1 1 leaf or canopy nitrogen concentration and light use efficiency at the ecosystem scale (Kergoat et 

12 al. 2008, Ollinger et al. 2008), However, we believe that remote sensing techniques that evaluate 

13 chlorophyll rather than nitrogen could have even greater potential for estimating ecosystem light 

14 use efficiency and GPP. 

15 From a biochemical perspective, only the PAR absorbed by photosynthetic pigments 

16 (designated as APAR c hi) enables photosynthetic processes, whereas the PAR absorbed by non- 

17 photosynthetic components such as boles, branches, stems, and litter is not used for CO 2 fixation. 

18 We designate chlorophyll-based FAPAR here as FAPAR chi . By definition, APAR canopy (the 

19 product of FAPAR canop y and PAR) is greater than APAR ch i (the product of FAPAR C hi and PAR). 

20 For linking to remote sensing applications, estimates of APAR ch i should provide more realistic 

21 GEP to wer and LUE tower values than similar estimates using APAR canop) ,. We define LUE based on 

22 APAR C hi versus APAR ca nopy as follows: 

23 
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LUE chl = 


GEP„ 


GEP 


FAPAR (hl * PAR APAR chl 


(3) 


LUE, 


GEP, 


GEP 


F APAR canopy * PAR APAR canopy 


(4) 


3 In earlier studies (Zhang et al., 2005, 2006), an approach to estimate FAPAR ch i was 

4 proposed using daily MODIS data. Since then, we have refined our algorithm to retrieve 

5 FAPAR C hi from MODIS imagery using the modified PROSPECT-SAIL2 model, PROSAIL-2 

6 (Zhang et al. 2005, 2006). The new version of this algorithm provides a statistical distribution of 

7 likely FAPARchi values for each cloud-free MODIS observation. 

8 In this article, we combine five years of flux, meteorological, and remote sensing data 

9 from a boreal aspen flux site to attain the following four objectives: (1) to present a method for 

10 estimating FAPAR chl and FAPAR canopy using single-date, scaled-up MODIS observations; (2) to 

1 1 apply the FAPAR ch i and FAPAR can opy algorithms to MODIS data acquired for 2001-2005 over 

12 this aspen flux site in Saskatchewan; (3) to link our estimates of MODIS FAPAR C hi and 

13 FAPARcanopy to the tower-based observations of PAR and GEP so as to derive LUE on both a 

14 unit cholorophyll area basis (LUE ch i, Eq.3 above) and for the whole canopy (LUE tanopy , Eq.4 

15 above); and (4) to compare our LUE ch i, LUE canopy and tow'er-based LUE estimates (i.e., LUE tower ) 

16 to see if the LUE C hi could provide a more physiologically realistic input to land surface process 

17 models. For this latter objective, we test the hypotheses that: (i)LUE canopy = LUE ch i; (ii) LUE ch i 

1 8 = LUEtower? and (iii) LUEcanopy = LUEtower- 

19 
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2 . 


Data and site descriptions 


2 2.1 Southern Old Aspen 

3 The Southern Old Aspen forest (SOA) was established in 1919 after a forest fire in Prince 

4 Albert National Park at the southern edge of the Canadian boreal forest (Barr et al. 2007). The 

5 eddy flux tower site (53.7°N, 106.2°W, 600 m elevation) is located ~50 km northwest of Prince 

6 Albert, Saskatchewan. SOA originated as part of the BOReal Ecosystem Atmosphere Study 

7 (BOREAS), and has continued operations under the Boreal Ecosystem Research and Monitoring 

8 Sites (BERMS) project and the Canadian Carbon Program, formerly the Fluxnet-Canada 

9 Research Network (FCRN). The vegetation around the tower site is primarily deciduous forest 

10 dominated 90% by aspen ( Populus tremuloides Michx.) with 10% balsam poplar ( Populus 

1 1 balsamifera L.) and a dense understory approximately 2 m tall of hazelnut ( Corylus cornuta 

12 Marsh) interspersed with green alder ( Alnus crispa Pursch). The aspen stand extends for at least 

13 3 km in all directions from the tower. The soil is mainly an orthic grey luvisol with a ~9 cm 

14 surface organic layer. The terrain is basically level. Mid-growing season leaf area index (LAI) 

15 under the tower varied from ~3.5 to 5.5. The climate is warm in summer (av. 17.5°C in July) and 

16 cold in winter (av. -19.1°C in January). Average annual precipitation was 412 mm during 2001- 

17 2005, whereas the long-term 1951-1980 average annual precipitation was higher at 484 mm 

18 (Griffis et al. 2003). A detailed analysis of the inter-annual variability of climatic factors at this 

19 site is presented in Barr et al. (2007). 

20 

21 2.2 GEPData 

22 GEP at SOA was estimated as the sum of Net Ecosystem Production (NEP) and 

23 ecosystem respiration. Respiration was modeled from night time and cold season NEP using soil 
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1 temperature. Details on CO 2 flux measurement methodology and obtaining estimates of GEP 

2 including gap filling can be found in Krishnan et al. (2006) and Barr et al. (2007). 

3 

4 2.3 Daily MODIS data 

5 Four MODIS daily products (v004) were used in this study: [1 ] surface reflectance 

6 (MOD09GHK and MYD09GHK); [2] observation viewing and illumination geometry 

7 (MODMGGAD and MYDMGGAD); [3] observation pointer information (MODPTHKM and 

8 MYDPTHKM); and [4] reflectance data quality descriptors (MOD09GST and MYD09GST). 

9 The MODIS imagery is nominally acquired with 500 m x 500 m spatial resolution at nadir, and 

10 >500 m spatial resolution for off-nadir views. The MODIS daily land surface reflectance 

1 1 product provides seven (of 36) spectral bands: red (620-670 nm, band 1), blue (459-479 nm, 

12 band 3), green (545-565 nm, band 4), near infrared (NIRi, 841—875 nm, band 2; NIR 2 , 1230— 

13 1250 nm, band 5), and short-wave infrared (SWIR|, 1628-1652 nm, band 6; SWIR 2 , 2105-2155 

14 nm, band 7). 

15 The MODIS daily observation viewing geometry product provides viewing and 

16 illumination geometry information (view zenith angle, VZA; view azimuth angle; sun zenith 

17 angle; and sun azimuth angle) at a nominal 1-km scale. The MODIS daily observation pointer 

1 8 product provides a reference, at a nominal 500 m scale, linking observations that intersect each 

19 pixel in the daily surface reflectance product to those given in the daily observation viewing 

20 geometry product. The MODIS daily reflectance data quality product provides summary quality 

21 information about MODIS daily surface reflectance conditions, including clouds, cloud shadow, 

22 land and water designations, aerosols, fire, snow, ice and bidirectional reflectance distribution 

23 function (BRDF) corrections, etc. All the MODIS data products are freely available at USGS 
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1 Earth Observing System Data Gateway (htt p://edcimswww.cr.usgs.gov/pub/imswelcome/) . and 

2 are delivered to users in a tile fashion, where each tile covers an area of 10° (latitude) by 10° 

3 (longitude). The software developed by the MODIS land team (MODLAND Tile Calculator 

4 http://modland.nascom.nasa.gov/cgi-bin/developer/tilemap.cgi) was utilized to determine the 

5 location of the SOA tower site in the MODIS products, including tile, row and column numbers. 

6 

7 3. Methods 

8 3.1 Spatial integration of MODIS data 

9 A possible source of discrepancy between tower-based photosynthesis and the MODIS 

10 standard GPP is the way the MODIS reflectance products are used to calculate FAPAR canopy 

1 1 and the MODIS standard GPP (Justice et al. 1998; Wolfe et al. 1998). This occurs because the 

12 areal coverage of the MODIS products used in the GPP calculations are not constant over the 

1 3 growing season and may not match the footprint of the flux tower site. The MODIS observations 

14 made at multiple times over a target actually cover somewhat different ground areas due to shifts 

15 in the ground track of the satellite, but are gridded into single, fixed grid cells. Footprints for off- 

16 nadir observations are increasingly larger and oblong in shape as VZAs increase. For example, 

17 areas associated with ground targets for imagery acquired with VZAs greater than 65° are at least 

18 nine times larger than those viewed by nadir observations (Wolfe et al. 1998). Consequently, 

19 both the standard MODIS FAPAR canop y and GPP products represent somewhat different, though 

20 overlapping or adjacent, areas when viewed frequently over time with different geometries. 

21 For the current analysis, we developed a method to scale-up the MODIS land band 

22 observations to 2500 m x 2500 m regions. We assumed that the scaled-up satellite 2500 m data 

23 and tower flux-based data should follow a similar pattern, although they are not necessarily 
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1 identical. This is a reasonable assumption since the flux footprint for the site described in the 

2 current study extends for at least three km in all directions from the flux tower. 

3 We acquired daily MODIS data (tile HI 1 V04) for 2001-2005. An example is given 

4 (Figure 1) for MODIS daily NIRi reflectances across the 5x5 grid area (where each grid cell is 

5 nominally 500 m) for nadir data on day of year (DOY) 224 in 2001. The MODIS relative 

6 reflectance data quality descriptor indicates that the reflectances were of high quality. 

7 We used 5x5 scaled up MODIS observations to produce similar ground sectors within the 

8 2.5 km x 2.5 km block area of the aspen forest (as in Fig. 1). A time series for the 5x5 block was 

9 created from MODIS daily data using the following criteria: (i) only observations that fell within 

10 the block were selected; (ii) an observation was excluded from consideration if the reflectance 

1 1 quality product indicated any quality problem; (iii) observations were averaged only if their 

12 geometries for view and illumination angles differed by less than five degrees; (iv) observations 

13 from different swaths were not mixed; and (v) the inclusion of at least 18 of the 25 grid cells 

14 were required to produce a scaled-up average observation for use in the subsequent analysis 

15 (Table I). 

16 The NDVI time series (Tucker 1979, Eq. 5 below) over five years (2001-2005) were used 

17 to determine average core growing season length where: 

18 

Pn!R. ~ P red ( C\ 

19 NDVI = 1 . (5) 

PnIR x T P red 

20 Other remote sensing indices were explored, e.g. the enhanced vegetation index (EVI; Huete et 

21 al. 1994) and the land surface w'ater index (LSWI, Xiao et al. 2005), but NDVI gave more 

22 coherent seasonal results (data not shown). Degree 6 polynomials were used to fit average 

23 seasonal curves to the five-year data collections for the MODIS NDVI and the mid-day tower 
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1 GEP, respectively. Derivative analyses were applied to the fitted curves. We used the first, 

2 second, and third derivatives to identify transition points in continuous data curves (Vina et al. 

3 2004) to retrieve the DOY in the five-year collections when seasonal changes occurred. We 

4 defined the beginning date of the growing season as the date when the second derivative reached 

5 its first local extreme value during the spring green-up period. We defined the end date for the 

6 growing season as that date when either the first or third derivatives attained a local extreme 

7 value at the end of summer (DOY > 255 at this site), initiating autumn senescence. Based on 

8 these analyses, we used NDVI for examining LUE and FAPAR dynamics for the 

9 photosynthetically active period, defined as occurring between DOY 152-259. We also defined 

10 the earlier and later dates for a second category representing periods of lower physiological 

1 1 activity (DOY between 121 and 151 and DOY between 260 and 287). 

12 

13 3.2 Estimating vegetation canopy characteristics 

14 3.2.1 Description of the PROSAIL-2 model 

15 The canopy- leaf-stcm-background coupled radiative transfer model, PROSAIL-2, is an 

16 updated version of the model used in earlier studies (Zhang et al. 2005, 2006). PROSAIL-2 

17 resulted from the combination of the PROSPECT model used to describe leaf characteristics in 

1 8 the canopy model, and the radiative transfer model, SAIL2. In the current version of SAIL2, we 

19 revised the expression for background and stem characteristics. Equations 5 & 6 in Zhang et al. 

20 (2006) were used to simulate soil and stem reflectance. In the present study, we used in situ 

21 measurements from the Oak Ridge National Laboratory website: 

22 (htt p://www.daac.oml.gov/BOREAS/boreas home page.html) to provide the “search ranges” for 

23 background (referred to as “back”) and stem spectral reflectance (referred to as “stem”). 
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1 Pback (^) = Pback . min U) + BACK A • (P bact ^ (A) ~ p, Mckmm (A)) (6) 

2 Pstcn, U) = Amm.min 0*) + A * M “ Ate*.™. &)) ’ ( 7 ) 

3 where 1 is the spectral wavelength, p is reflectance, BACKa and STEMa are variables 

4 describing reflectance values for background and stem. We used the maximum and minimum 

5 reflectance values of background and stem as their upper and lower value limits. The fourteen 

6 free variables used in the PROSAIL-2 for this study are summarized in Table 2. Five variables 

7 were used to describe leaf characteristics: a leaf internal structure variable (N), leaf total 

8 photosynthetic pigment content (C ab ), leaf dry matter content (C m ), leaf water thickness (C w ), and 

9 leaf brown pigment content (C br own)- The top of canopy reflectance was composed of leaf, stem 

10 and background contributions. Stem fraction (SFRAC) and cover fraction (CF) were used to 

1 1 decompose leaf, stem and background components. Refer to Zhang et al. (2005, 2006) for more 

12 details on PROSAIL-2. 

13 

14 3.2.2 Description of the FAPAR c m and FAPAR canopy algorithm 

15 The variables in Tab. 2 have to be estimated to calculate FAPAR c hi and FAPAR canop y 

16 using PROSAIL-2. The Metropolis algorithm, a Markov Chain Monte Carlo (MCMC) method, 

17 was adopted to find solutions expressed as posterior distributions of the variables. The posterior 

18 distributions of the variables are the product of their prior distributions and the likelihood 

19 calculated with a Bayesian analysis (http://en.wikipedia.org/wiki/Posterior_probability). The 

20 prior distributions of the input variables were assumed to be uniform ( http://en.wikipedia.org/ 

21 wiki/Uniform_distribution_(continuous)). Each MODIS reflectance observation [p obs ] for the 

22 seven land bands (red, NIR], blue, green, NIR 2 , SWIRj and SWIRj), and associated VZA [0 V , in 

23 degrees], relative view azimuth angle [<j>, in degrees], and solar zenith angle [0 S , in degrees] 
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1 contains some noise, although small differences in angles may be ignored. We treat each SOA 

2 reflectance observation as a sample of the following distribution: 

3 

p ~ {p nhs (A, e v (1 + 3N(0,1)), e, (1 + 3N(0,1)), <p{] + 3N(0,1))) } • d + 0.05N(0,f)) (8) 

5 

6 where N(0,1) is the normal distribution with a mean of zero and SD = 1. 

7 We may use as many samples from the distribution (from equation 8) as we desire. For 

8 five of the spectral bands (red, green, NIRi, NIR 2 and SWIRi), we calculated the log-likelihood 

9 [using equations 1& 2 from Zhang et al. 2005] and then performed an acceptance test [using 

10 equation 3 from Zhang et al. 2005]. A new randomly generated “proposed” value was accepted 

1 1 only if it passed acceptance tests conducted on all five bands. The same adaptive algorithm 

1 2 [using equation 4 from Zhang et al. 2005] was used to accelerate the speed of convergence of the 

13 MCMC algorithm. With the posterior distributions of the variables (Tab. 2) now calculated, we 

14 forward-simulated the fractions of APAR for canopy, leaf, and photosynthetic pigments (Zhang 

15 et al. 2005, 2006) for each MODIS 5x5 scaled-up observation that met our quality rules during 

16 the five-year period, using the PROSAIL-2 model. The product of LAI, CF, and leaf 

17 photosynthetic pigment (jlg/cnT) describes the average photosynthetic pigment content for each 

18 5x5 MODIS aspen forest observation. Similarly, the product of LAI, CF and leaf water content 

19 (cm) describes the average water content of vegetation and, the product of LAI, CF and leaf dry 

20 matter (g/cm 2 ) provides the average dry matter content of vegetation in each MODIS observation. 

21 
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1 3.2.3 Calculation of AP ARchb APARcanopy^ ElJEcht EEJEcanopy and ELJEf^^cr 

2 The PAR absorbed by the whole SOA canopy (APAR ca nopy) was determined as the 

3 product of PAR and the median value of the MODIS-derived FAPAR can opy distribution (see 

4 section above). Likewise, APAR ch i for only the photosynthetic pigments of the foliage 

5 component was determined as the product of PAR and the median value of the MODIS-derived 

6 FAPARchi distribution. For these calculations, we used the GEP and average incident canopy 

7 photosynthetic photon flux density (PPFD) measured with Li- 1 90S A PAR sensors (Licor Inc., 

8 Lincoln, NE) over the SOA tower site at 90-minute intervals centered on the satellite overpass 

9 time. LUE for the foliage component of the forest at SOA (LUE ch i, Eq.3) was computed as the 

10 ratio of GEP to APAR c hi. LUE for the whole SOA forest canopy (LUE canopy , Eq.4) was 

1 1 computed as the ratio of GEP to APAR can0 p y . The LUE tower (Eq. 2 above) data used in this paper 

1 2 are cited from Krishnan et al. (2006) where they defined LUE t0W er as the ratio 

13 GEP t0W er/APAR t0W er where APARtower was estimated using Eq. (1) in Barr et al. (2007) that was 

14 based on measured downwelling and upwelling PAR, overstory and understory clumping indices, 

15 measured stem area indices, and the estimated daily LAI. Therefore, LUE can opy and LUE ch i use 

16 MODIS-derived FAPARs whereas LUE towe r derives its FAPAR equivalent from site-based 

17 radiometer and structural measurements. The GEP and PAR variables, on the other hand, are the 

18 same for the three LUE calculations. 

19 The REGRESS function in MATLAB was used to statistically analyze regression 

20 relationships between: (i) GEP and APAR ch i; (ii) GEP and APAR canopy ; (iii) annual average 

2 1 LUE C hi ; and annual average LUE canopy over the five-year core growing season period: and (iv) 

22 annual average LUE towe r and annual average LUE ch i. We also used the student t-test function in 

23 MATLAB to test if the slope between GEP and APAR C hi and the slope between GEP and 
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1 APARcanopy were significantly different. The same t-test function was also applied to test if the 

2 annual average LUE C hi- LUE can opy and LUE t0Wcr time series over the entire five-year period were 

3 significantly different. 

4 

5 3.2.4 Controls on inter-annual LUE C hi and LUE C anopy 

6 To determine the controls on the inter-annual variation of LUE, we analyzed the 

7 relationship between LUE c hi, LUE canopy , soil w'ater content (SWC), and precipitation as well as 

8 the three factors that determine canopy chlorophyll concentration, i.e., average leaf chlorophyll 

9 concentration, leaf area index (LAI), and cover fraction (CF). Soil water content was measured 

10 with eight time domain reflectometry (TDR) probes (Moisture Point type B, Gabel Corp., 

1 1 Victoria, Canada) placed at 10 m intervals with measurements at depths of 0-15, 15-30, 30-60, 

12 60-90, and 90-120 cm, although SWC values were used only for the 0-30 cm zone because this is 

13 where 90% of the roots are found (Barr et al. 2007). Rain precipitation was measured with a 

14 Geonor T200 accumulation rain gauge (Geonor Inc, Milford, PA) supplemented by a CS700 

15 tipping bucket rain gauge (Campbell Scientific Inc, Edmonton, AB). The half-hourly 

16 precipitation measurements at the SOA tower site in 2001-2005 were summed into monthly rain 

17 precipitation totals, from which annual values were 235 mm, 285.8 mm, 261 mm, 667 mm and 

18 614 mm for 2001 - 2005, respectively. 

19 

20 4. Results 

21 4.1 PROSAIL-2 Derived Canopy Variables including APAR C hi and APAR canopv 

22 The posterior distributions of SOA canopy variables from the PROSAIL-2 model display 

23 seasonal variation, as shown for nine of the seventeen MODIS daily 5x5 observations in 2005 
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1 (Table 3). The other observations had similar distributions as those shown in Tab. 3 so they are 

2 not presented. Several canopy variables are shown: LAI, CF, total photosynthetic pigment 

3 content, water content, dry matter content, FAPAR canop y, and FAPAR c hi. Differences in total 

4 chlorophyll concentration were related to changes in all three of the factors used in its calculation, 

5 i.e., average leaf chlorophyll concentration, LAI, and CF. The considerably higher FAPAR can opy 

6 (0.47 - 0.87), as compared to FAPAR chl (0.03 - 0.70), results from PAR absorption primarily by 

7 non-photosynthetic canopy components. The average APAR ch i at SOA over the five-year period 

8 was roughly 65% of APAR can opy (slope of the all-data relationship in Fig. 2). We also found that 

9 the average ratios for APAR C ht : APAR canopy were different between the core growing season, 

10 DOY 152 - 259 (♦, slope = 0.71); and the combined early and late periods of the season, DOY 

1 1 <152 and DOY>259 (A, slope = 0.34) (Fig. 2). The correlation between APAR ch i and 

12 APAR can opy showed the highest correlation (r = 0.87) during the core growing season, DOY 152- 

13 259 (Fig. 2). 

14 

15 4.2 LUE C hi and LUE canony Over the Five-Year Period 

16 Average values for LUE can0 py and LUE c hi over the five-year period (2001-2005) were 

17 0.0155 and 0.0241 (imol C ptnof 1 PPFD, respectively, as determined from the slopes of the 

18 GEP:APAR relationships (Figure 3 for APAR canopy ; Figure 4 for APAR ch i, p<0.0001). For the 

19 entire study period (DOY ranging from 121 to 287), there was a stronger correlation between 

20 GEP and APARchi, (Fig. 4, r 2 = 0.78) compared GEP:APAR can0 p y (Fig. 3, r 2 = 0.33). The 95% 

21 confidence intervals for the five-year average LUE canopv and LUE cW did not overlap, i.e., they 

22 ranged from 0.0141 to 0.0169 and from 0.0229 to 0.0253, respectively. 
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1 For the core growing season only (DOY 152 - 259), the five-year average values of 

2 LUE can0 py and LUE ch i were 0.0173 and 0.0243 |iimol C pmol" 1 PPFD, respectively (Figure 3 for 

3 APARcanopy* Figure 4 for APAR C hi, p<0.0006). Once again, the 95% confidence intervals for the 

4 five-year average LUE canop y and LUE C hi did not overlap, i.e., they ranged from 0.0163 to 0.0183 

5 and from 0.0230 to 0.0256, respectively. 

6 APAR c hi and LUE C hi captured more seasonal variation than their whole canopy 

7 counterparts (APAR canopy , LUE can0 py) (Figs. 3 and 4). For example, during the early and late 

8 growing season, average LUE c hi over five years was 0.0208 pmol C pmol 1 PPFD which was 

9 lower than the core growing season value of 0.0243 pmol C pmol' 1 PPFD (Fig. 4). 

10 

11 4.3 Inter-Annual Variability in LUE 

12 There were also inter-annual variations for APAR c hi and APAR canopy (mean ± SE, the 

13 standard error) with consistently and significantly lower values for APAR ch i compared to 

14 APAR can opy during the core growing season (Figure 5) (p<0.005). Consequently, LUE ch i was 

15 consistently and significantly higher than LUE can opy (Figure 6) (p<0.005) and LUE t hi in 2004- 

16 2005 was higher than for the three earlier years during the growing season. The effect of the 

1 7 2003 drought was apparent, such that the average LUE t hi and LUEcanopy in 2004 were much 

18 higher than in 2003. The annual means of core growing season LUE c hi (pmol C pmoF 1 PPFD) 

19 (± SE) were: 0.0242 ± 0.0012 (2001), 0.0245 ± 0.0015 (2002); 0.0225 ± 0.0018 (2003); 0.0310 ± 

20 0.0022 (2004); and 0.0267 ± 0.0019 (2005) (Tab. 4). The maximum LUE (£ max ) value for 

2 1 broadleaf deciduous forests set by the biome look-up table of the MOD1S photosynthesis model 

22 is 0.0203 pmol C pmol' 1 PPFD (Heinsch et al., 2003). Four of the five annual LUE canopy values 
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that we calculated were lower than e max . However, the 2004 LUE canopy was much higher than the 


2 MODIS £ max value (Fig. 6). 

3 Average LAI, CF, canopy chlorophyll concentration, FAPAR canopy , FAPAR ch i, LUE canop y 

4 and LUE C hi during core growing season (152< DOY< 259) for each individual year are shown in 

5 Table 4. Inter-annual variations of canopy chlorophyll concentration were influenced by all the 

6 three factors (r - = 0.91 for leaf chlorophyll concentration; r = 0.32 for LAI; and r = 0.49 for CF). 

7 Average annual canopy chlorophyll concentrations were correlated with both FAPAR ch i (r 2 = 

8 0.8 1 ) and APAR ch i (r 2 = 0.63). 

9 

10 4.4. Comparisons with LUE, ower 

1 1 Whereas LUE canop y differed significantly from LUE t0W cr for each of the five years 

12 (p<0.001), LUEchi was essentially the same as LUEtower (pLO.47, Figure 6). 

13 

14 5, Discussion 

15 LUEchi captured more seasonal (Tab. 3) and inter-annual variation than LUE can0 p y and 

16 provided an improved overall relationship to GEP (Figs. 3 and 4). Krishnan et al. (2006) reported 

17 that the annual average LUE t0W er at the SOA tower for 2001-2005 was 0.0229 - 0.0302 pmol C 

18 pmol' 1 PPFD. Their in situ average LUE toW er estimate matched well with our average MODIS- 

19 derived LUE c hi (0.0229 - 0.0302 vs. 0.0225 - 0.0310 pmol C pmol ' PPFD from Fig. 4) but not 

20 with the MODIS-derived LUE can opy over five growing seasons. Furthermore, the in situ tower- 

2 1 based LUE t0 wer and the MODIS-derived LUE C hi were also similar in each of the five years, 

22 exhibiting their highest values in 2004-2005 (Fig. 6). 
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1 The annual LUE ca no P y values substantially underestimated the tower-based estimates 

2 (Krishnan et al. 2006) (Fig. 6). It is interesting to note that the five-year average LUE catl0 p y for 

3 the mid-growing season period was higher than LUE canop y for the early and late season periods 

4 (Fig. 3). In comparison, the five-year average LUE C hi for all data was consistent with LUE C hi for 

5 the mid-growing season period, and also close to LUE ch i for the early and late season periods 

6 (Fig. 4). Additionally, the maximum LUE value (£ max ) from the MODIS biome look-up table for 

7 broadleaf deciduous forests (Heinsch et al. 2003) tended to be significantly lower than both 

8 LUE C hi and LUE t0 wer (Fig. 6). This study has demonstrated that (1) LUE ch i values are more 

9 comparable to ground-based observations of LUE toW er than LUEcanopy! (2) APAR C hi values are 

10 more comparable to ground-based observations of APAR tower when the sky is clear than 

11 APAR ca nopy; (3) FAPAR c hi values derived from MODIS observations are more realistic and 

12 useful for estimating of APAR lower when the sky is clear than FAPAR canopy values. The 

13 conclusions are supported by measurements and simulations of photosynthetic vs. non- 

14 photosynthetic vegetation by Chen et al. (2006) at other flux sites in the region. Consequently, 

15 we accept the hypothesis that LUE C hi = LUE toW er, whereas we reject the other two hypotheses that 

1 6 LUEcanopy = LUEchi Or that LUEcanopy = LUE| 0wer . 

17 A prolonged three-year drought began in late summer of 2001 and the most severe 

18 conditions occurred in 2003 (Barr et al. 2007), which also produced the lowest LUE c hi (Fig. 6). 

19 There were eleven consecutive months prior to August 2003 when precipitation was <50 mm per 

20 month. The water table depth decreased from 3 m in 2001 to 4 m in 2003 (Barr et al. 2007). Low 

21 precipitation in 2001 through 2003 caused soil water content in the shallow layer (0-0.15 m) to 

22 begin to drop in August 2001 and it kept dropping through 2002 and 2003. The lowest soil water 

23 content occurred in 2003 and was one-third lower than the pre-drought mean value and was close 
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1 to the permanent wilting point. Surface conductance declined during the drought years and 

2 reached its lowest value in 2003 (Krishnan et al. 2006). Figure 7 compares the average LUE C hi 

3 and LUE can0 py during DOY 152-259 of each year (2001 - 2005) with the cumulative rain 

4 precipitation during that time period. This analysis suggests that the cumulative precipitation 

5 during 2001-2005 may have had a significant influence on LUE. Increased rain precipitation in 

6 2004 and 2005 recharged the soil and increased LUE and our MODIS-based LUE ch i estimates 

7 were sensitive to this phenomenon. 

8 FAPARchi and LUE ch i are more physiologically realistic ways of quantifying the PAR 

9 absorbed and used for photosynthesis. These are integrative measures that inherently account for 

10 some of the impacts of mid- to long-term stresses since they also reflect changes in leaf area and 

1 1 cover fraction. For example, inter-annual variability of total canopy chlorophyll concentration 

12 was influenced by changes in average leaf chlorophyll concentration, LAI and CF (Table 4). 

13 The 2003 drought lowered all three of these factors at our study area relative to 2001 and 2002. 

14 Thus, chlorophyll-based measures of APAR and LUE have the potential to more directly account 

15 for environmental limitations and thus reduce the impact of the uncertainty in estimating 

16 temperature, vapor pressure deficit, and soil moisture for specific pixels. Waring et al. (1995) 

17 found a strong correlation between upper canopy leaf chlorophyll concentration of the major 

1 8 hardwood species and maximum light use efficiency at the primarily deciduous Harvard Forest 

19 flux site. Their finding supports the idea of a direct link between pigment concentration and 

20 LUEtower for deciduous forests. Ollinger et al. (2008) and Kergoat et al. (2008) showed a 

21 significant positive relationship between whole canopy nitrogen concentration and canopy 

22 maximum LUE but the relationship between canopy nitrogen and chlorophyll was not described 

23 and they did not examine inter-annual or seasonal variability at any of their sites. Their results 
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1 lead us to believe that an inter-site analysis based on APAR C hi and LUE C hi using our methodology 

2 would yield even stronger relationships. 

3 Chlorophyll-based measures of APAR and LUE, however, will not account for 

4 limitations due to short-term environmental extremes so modulation of light use efficiency by 

5 environmental stresses will still need to be considered. Furthermore, photosynthesis of evergreen 

6 conifer forests are less sensitive to changes in chlorophyll than are broadleaf forests as evidenced 

7 by the continuous green color of conifer forests even during the coldest periods of winter. As 

8 well, forests that have attained maximum height are subject to significant hydraulic and stomatal 

9 limitations (Ryan et al. 2004, 1997) that may or may not be reflected in their total canopy 

10 chlorophyll concentration. Brodribb and Feild (2000) demonstrated a highly significant 

1 1 correlation between hydraulic conductivity, maximum photosynthetic capacity, and quantum 

12 yield 23 rain forest species. 

1 3 To our knowledge, this study represents the first time that LUE C hi has been estimated by 

14 linking tower flux data, a biophysical radiative transfer model, and satellite spectral observations. 

15 Our predicted values for LUE c hi agree well with in situ data in respect to both amplitude and 

16 seasonal phase; and our modeled LUE ch i successfully described the actual dynamics captured by 

17 tower fluxes. We believe that it could be useful to couple this type of FAPAR C hi approach into 

18 regional/global carbon cycle models, land surface process models and general circulation models. 

19 This paper has demonstrated some of the possible benefits of using FAPAR C hi as an 

20 operational data product for carbon cycle modeling. The hyperspectral satellite sensors that are 

21 in orbit (e.g., EO-1 /Hyperion) or currently under development (e.g., HyspIRI) could help us 

22 obtain even more robust estimates of FAPAR c hi due to the greater sensitivity of these sensors to 

23 pigment concentrations and other biochemical properties of foliage (Coops et al. 2002). A 


23 



1 narrower field of view 7 could also be helpful for resolving chlorophyll dynamics at scales more 

2 representative (e.g., <100 m) of the spatial structure typical of most forest stands. 

3 

4 6. Acknowledgments 

5 This project was supported by the NASA Carbon Cycle Science Program (Dr. Diane Wickland, 

6 Program Manager) under the auspices of the North American Carbon Program. Support was 

7 provided to the FCRN and BERMS programs through several Canadian funding sources 

8 (Environment Canada, the Canadian Forest Service, NSERC, CFCAS, BIOCAP, Action Plan 

9 2000 and PERD). 

10 

1 1 7. References 

12 Barr, A.G., Black, T.A., Flogg, E.H., Griffis, T.J., Morgenstem, K., Kljun, N., Theede, A., & 

13 Nesic, Z. (2007). Climatic controls on the carbon and water balances of a boreal aspen 

14 forest, 1994-2003. Global Change Biology, 13, 561-576. 

15 Brodribb, T.J. & Feild, T.S. (2000). Stem hydraulic supply is linked to leaf photosynthetic 

16 capacity: evidence from New Caledonian and Tasmanian rainforests. Plant, Cell and 

17 Environment 23, 1381-1388 

18 Coops, N.C., Smith, M.-L., Martin, M.E., Ollinger, S.V., & Held, A. (2002). Predicting eucalypt 

19 biochemistry from HYPERION and HYMAP imagery. IGARSS'02, 2, 790 - 792 

20 Gower, S.T., Kucharik, C.J., & Norman, J.M. (1999). Direct and indirect estimation of leaf area 

2 1 index, fAPAR and net primary production of terrestrial ecosystems. Remote Sensing of 

22 Environment, 70, 29-5 1 


24 



2 


Griffis, T.J., Black, T.A., Morgenstem, K„ Barr, A.G., Nesic, Z„ Drewitt, G.B., Gaumont-Guay, 
D., & McCaughey, J.H. (2003). Ecophysiological controls on the carbon balances of 

3 three southern boreal forests. Agricultural and Forest Meteorology, 117 , 53-71 

4 Heinsch, F.A., M. Reeves, P. Votava, S. Kang, C. Milesi, M. Zhao, J. Glassy, W.M. Jolly, R. 

5 Loehman, C.F. Bowker, J. Kimball, R. Nemani, S.W. Running. (2003). User's Guide: 

6 GPP and NPP (MOD 17A2/A3) Products. NASA MODIS hand Algorithm. Version 2.0. 

7 NASA Algorithm Theoretical Basis Document (ATBD). 

8 Huete, A., Justice, C., & Liu, H. (1994). Development of Vegetation and Soil Indexes for Modis- 

9 Eos. Remote Sensing of Environment, 49, 224-234 

10 Justice, C.O., Vermote, E., Townshend, J.R.G., Defries, R., Roy, D.P., Hall, D.K., Salomonson, 

1 1 V.V., Privette, J.L., Riggs, G., Strahler, A., Lucht, W., Myneni, R.B., Rnyazikhin, Y., 

12 Running, S.W., Nemani, R.R., Wan, Z.M., Huete, A.R., van Leeuwen, W., Wolfe, R.E., 

13 Giglio, L., Muller, J.P., Lewis, P., & Barnsley, M.J. (1998). The Moderate Resolution 

14 Imaging Spectroradiometer (MODIS): Land remote sensing for global change research. 

15 IEEE Transactions on Geoscience and Remote Sensing, 36, 1228-1249 

16 Kergoat, L., Lafont, S., Ameth, A., Le Dantec, V., & Saugier, B. (2008). Nitrogen controls 

17 plant canopy light-use efficiency in temperate and boreal ecosystems. Journal of 

1 8 Geophysical Research, 113, G0401 7. doi : 1 0. 1 029/2007 JG000676. 

19 Krishnan, P., Black, T.A., Grant, N.J., Barr, A.G., Hogg, E.H., Jassal, R.S., & Morgenstem, K. 

20 (2006). Impact of changing soil moisture distribution on net ecosystem productivity of a 

2 1 boreal aspen forest during and following drought. Agricultural and Forest Meteorology, 

22 1 39, 208-233 


25 



1 Landsberg, J.J., & Waring, R.H. (1997). A generalized model of forest productivity using 

2 simplified concepts of radiation-use efficiency, carbon balance and partitioning. Forest 

3 Ecology and Management, 95, 209-228 

4 Law, B.E., Waring, R.H., Anthoni, P.M., & Aber, J.D. (2000). Measurements of gross and net 

5 ecosystem productivity and water vapour exchange of a Pinus ponderosa ecosystem, and 

6 an evaluation of two generalized models. Global Change Biology, 6 , 155-168 

7 Medlyn, B.E. (1998). Physiological basis of the light use efficiency model. Tree Physiology, 18, 

8 167-176 

9 Mohr, H., Schopfer, P., Lawlor, G., & Lawlor, D.W. (Eds.) (1995). Plant Physiology. Springer 

10 Monteith, J.L. (1972). Solar-Radiation and productivity in tropical ecosystems. Journal of 

1 1 Applied Ecology, 9, 747-766 

12 Monteith, J.L. (1977). Climate and efficiency of crop production in Britain. Philosophical 

13 Transaction of the Royal Society of London B: Biological Sciences, 281, 271-294 

14 Ollinger, S.V., Richardson, A.D., Martin, M.E., Hollinger, D.Y., Frolking, S.E., Reich, P.B., 

15 Plourde, L.C., Katul, G.G., Munger, J.W., Oren, R.O., Smith, M.-L., Paw U, K.T., 

16 Bolstad, P.V., Cook, B.D., Day, M.C., Martin, T.A., Monson, R.K., & Schmid, H.P. 

17 (2008). Canopy nitrogen, carbon assimilation, and albedo in temperate and boreal forests 

18 Functional relations and potential climate feedbacks. Proceedings, National Academy of 

19 Sciences, 105, 19336-19341. 

20 Prince, S.D., Goetz, S.J., & Goward, S.N. (1995). Monitoring Primary Production from Earth 

21 Observing Satellites. Water Air and Soil Pollution, 82, 509-522 

22 Prince, S.D., & Goward, S.N. (1995). Global primary production: A remote sensing approach. 

23 Journal of Biogeography, 22, 815-835 


26 



Prince, S.D., & Goward, S.N. (1996). Evaluation of the NOAA/NASA Pathfinder AVHRR Land 


2 Data Set for global primary production modelling. International Journal of Remote 

3 Sensing, 17, 217-221 

4 Prince, S.D., Goward, S.N., Goetz, S„ & Czajkowski, K. (2000). Interannual Atmosphere- 

5 Biosphere Variation: Implications for observation and modeling. Journal of Geophysical 

6 Research-Atmospheres, 1 05, 20055-20063 

7 Ruimy, A., Dedieu, G., & Saugier, B. (1996a). TURC: A diagnostic model of continental gross 

8 primary productivity and net primary productivity. Global Biogeochemical Cycles, 10, 

9 269-285 

10 Ruimy, A., Kergoat, L., Field, C.B., & Saugier, B. (1996b). The use of C02 flux measurements 

1 1 in models of the global terrestrial carbon budget. Global Change Biology, 2, 287-296 

12 Ruimy, A., Saugier, B., & Dedieu, G. (1994). Methodology for estimation of terrestrial net 

13 primary production from remotely sensed data. Journal of Geophysical Research - 

14 Atmospheres, 99, 5263-5283 

15 Running, S., Nemani, R., Heinsch, F., Zhao, M., Reeves, M., & Hashimoto, H. (2004). A 

16 continuous satellite-derived measure of global terrestrial primary production. BioScience, 

17 54, 547-560 

18 Running, S.W., Baldocchi, D.D., Turner. D.P., Gower, S.T., Bakwin, P.S., & Hibbard, K.A. 

19 (1999a). A global terrestrial monitoring network integrating tower fluxes, flask sampling. 

20 ecosystem modeling and EOS satellite data. Remote Sensing of Environment, 70, 108-127 

21 Running, S.W., Justice, C.O., Salomonson, V., Hall, D., Barker, J., Kaufmann, Y.J., Strahler, 

22 A.H., Huete, A.R., Muller, J.P., Vanderbilt, V., Wan, Z.M., Teillet, P., & Cameggie, D. 


27 



1 (1994). Terrestrial Remote-Sensing Science and Algorithms Planned for Eos Modis. 

2 International Journal of Remote Sensing, 75, 3587-3620 

3 Running, S.W., Nemani. R„ Glassy, J.M., & Thornton, P.E. (1999b). MODIS daily 

4 photosynthesis (PSN) and annual net primary production (NPP) product (MOD17), 

5 Algorithm Theoretical Basis Document, Version 3.0 

6 Running, S.W., Thornton, P.E., Nemani, R., & Glassy, J.M. (2000). Global terrestrial gross and 

7 net primary productivity from the Earth Observing System. In O.E. Sala, R.B. Jackson, 

8 H.A. Mooney & R.W. Howarth (Eds.), Methods in Ecosystem Science (pp. 44-57). New 

9 York: Springer Verlag 

10 Ryan M.G., Binkley D., Fownes J.H. (1997) Age-related decline in forest productivity: pattern 

1 1 and process. Advances in Ecological Research, 27, 213-262. 

12 Ryan M.G., Binkley D„ Fownes J.H., Giardina C.P., Senock, R.S. (2004). An experimental test 

13 of the causes of forest growth decline with stand age. Ecological Monographs, 74, 393- 

14 414. 

15 Sellers, P.J., Los, S.O., Tucker, C .J., Justice, C.O., Dazlich, D.A., Collatz, G.J., & Randall, D.A. 

16 (1996a). A revised land surface parameterization (SiB2) for atmospheric GCMs .2. The 

17 generation of global fields of terrestrial biophysical parameters from satellite data. 

18 Journal of Climate, 9, 706-737 

19 Sellers, P.J., Randall, D.A., Collatz, G.J., Berry, J.A., Field, C.B., Dazlich, D.A., Zhang, C., 

20 Colleio, G.D., & Bounoua, L. (1996b). A revised land surface parameterization (SiB2) 

21 for atmospheric GCMs .1. Model formulation. Journal of Climate, 9, 676-705 

22 Tucker, C.J. (1979). Red and Photographic Infrared Linear Combinations for Monitoring 

23 Vegetation. Remote Sensing of Environment, 8, 127-150 


28 



1 Turner, D.P., Ollinger, S., Smith, M.L., Krankina, O., & Gregory, M. (2004). Scaling net 

2 primary production to a MOD IS footprint in support of Earth observing system product 

3 validation. International Journal of Remote Sensing, 25, 1961-1979 

4 Turner, D.P., Ritts, W.D., Cohen, W.B., Gower, S.T., Zhao, M.S., Running, S.W., Wofsy, S.C., 

5 Urbanski, S., Dunn, A.L., & Munger, J.W. (2003). Scaling Gross Primary Production 

6 (GPP) over boreal and deciduous forest landscapes in support of MODIS GPP product 

7 validation. Remote Sensing of Environment, 88, 256-270 

8 Turner, D.P., Ritts, W.D., Zhao, M., Kurc, S.A., Dunn, A.L., Wofsy, S., Small, E., & Running, 

9 S.W. (2006). Assessing interannual variation in MODIS-based estimates of gross primary 

10 production. Ieee Computational Science & Engineering, 44, 1899-1907 

1 1 Vina, A., A. Gitelson, D. Rundquist, G. Keydan, B. Leavitt, and J. Schepers (2004). Monitoring 

12 Maize (Zeo mays L.) Phenology with Remote Sensing. Agronomy Journal 96:1 139-1 147. 

13 Waring, R.H., Law, B.E., Goulden, M.L., Bassow, S.L., McCreight, R.W., Wofsy, S.C., Bazzaz, 

14 F.A. (1995). Scaling gross ecosystem production at Harvard Forest with remote sensing: 

15 a comparison of estimates from a constrained quantum-use efficiency model and eddy 

16 correlation. Plant, Cell, Environment 18, 1201-1213. 

17 Wolfe, R., Roy, D„ & Vermote, E. (1998). The MODIS land data storage, gridding and 

18 compositing methodology: Level 2 Grid. IEEE Transactions on Geosciences and Remote 

19 Sensing, 36, 1324-1338 

20 Xiao, X.M., Zhang, Q.Y., Hollinger, D., Aber, J., & Moore, B. (2005). Modeling gross primary 

2 1 production of an evergreen needleleaf forest using modis and climate data. Ecological 

22 Applications, 15, 954-969 


29 



1 


Yoder, B.J., & Waring, R.H. (1994). The Normalized Difference Vegetation Index of Small 
Douglas-Fir Canopies with Varying Chlorophyll Concentrations. Remote Sensing of 


2 

3 Environment, 49, 81-91 

4 Zhang, Q., Xiao, X., Brasw'ell, B.. Linder, E., Ollinger, S„ Smith, M.L., Jenkins, J.P., Baret, F., 

5 Richardson, A.D., Moore, B.. & Minocha, R. (2006). Characterization of seasonal 

6 variation of forest canopy in a temperate deciduous broadleaf forest, using daily MODIS 

7 data Remote Sens. Environ., 105, 189-203 

8 Zhang, Q.Y., Xiao, X.M., Braswell, B., Linder, E., Baret, F., & Moore, B. (2005). Estimating 

9 light absorption by chlorophyll, leaf and canopy in a deciduous broadleaf forest using 

10 MODIS data and a radiative transfer model. Remote Sens. Environ., 99, 357-371 

11 
12 

13 

14 


30 



List of abbreviations and acronyms 
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3 

BERMS - 

Boreal Ecosystem Research and Monitoring Sites 

4 

BOREAS - 

BOReal Ecosystem-Atmosphere Study 

5 

BRDF - 

bidirectional reflectance distribution function 

6 

chi 

chlorophyll 

7 

EVI 

enhanced vegetation index 

8 

FAPAR - 

fraction of absorbed photosynthetically active radiation 

9 

FCRN - 

Fluxnet-Canada Research Network 

10 

LSWI 

land surface water index 

11 

LUE 

light use efficiency 

12 

MCMC - 

Markov Chain Monte Carlo 

13 

MOD IS - 

MODerate resolution Imaging Spectroradiometer 

14 

NDVI 

normalized difference vegetation index 

15 

SOA 

the Southern Old Aspen tower site 


16 

17 

18 

19 

20 
21 
22 

23 

24 
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1 Table 1. Days of useable daily 5x5 MODIS obserx’ations over the Southern Old Aspen site (SOA) 

2 during 2001-2005 (n=76) 


Year 

Day of Year (DOY) 

2001 

144, 153, 155, 185, 189, 190, 215, 216, 218, 222, 

n=22 

224, 228, 233, 236, 245, 247, 249, 258, 261, 265, 


272, 278 

2002 

156, 174, 177, 179, 188, 190, 193, 197, 216, 232, 

n=16 

235,236, 239, 241,261,268 


139, 145, 148, 155, 166, 168, 197, 214, 226, 230, 


237, 246 

2004 

182, 199, 206, 207, 208, 212, 217, 226, 247 

n=9 


2005 

121, 122, 123, 128, 132, 150, 160, 173, 196, 219, 

n=17 

233, 244, 245, 247, 249, 285, 287 
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Table 2. A list of variables in the PROSAIL-2 model and the search ranges used for inversion. 



Variable 

Description 

Unit 

Search 

range 

Biophysical 

/biochemical 

PAI 

Plant area index, i.e., leaf +stem area 
index 

m 2 /m 2 

1-7.5 

variables 

SFRAC 

Stem fraction 


0-1 


CF 

Cover fraction: area of land covered by 
vegetation/total area of land 


0.5-1 


Cab 

Leaf photosynthetic pigments including 
chlorophyll a+b and carotenoids 

pg/cm 2 

0-150 


N 

Leaf structure variable: measure of the 
internal structure of the leaf 


1.0 -4.5 


c w 

Leaf equivalent water thickness 

cm 

0.001-0.15 


C m 

Leaf dry matter content 

g/cm 2 

0.001 - 0.04 


C brown 

Leaf brown pigment content 


0.00001 -8 


LFINC 

Mean leaf inclination angle 

degree 

10-89 


STINC 

Mean stem inclination angle 

degree 

10-89 


LFHOT 

Leaf BRDF variable: length of leaf/ 
height of vegetation canopy 


0 

1 

o 


STHOT 

Stem BRDF variable: length of stem/ 
height of vegetation canopy 


0-0.9 


STEMa 

Stem reflectance variable 


o 

i 

o 

o' 


BACKa 

Background reflectance variable 


p 

1 

o 

o 
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Table 3. Median values for photosynthetic pigments, water content, and dry matter at the grid 
cell level. LAI, CF, FAPAR ccmopy and FA PAR,/,/ were estimated using PROSAIL-2 from daily 5x5 
MODIS observations over the Southern Old Aspen site (SO A) in 2005. Note that the grid cell 
estimate = leaf level estimate*LAI*CF (where “estimate” refers to either photosynthetic 
pigments, water content, or dry matter, respectively ) 


DOY 

in 2005 

LAI 


Photosynthetic 

pigment 

(|4g/cm 2 ) 

Water 

content 

(cm) 

Dry 

matter 

(g/cm 2 ) 

F AP AR can0 py 
0-1 

FAPAR chl 

0-1 

132 

0.27 


3.75 

0.0005 

0.0004 

0.470 


150 

1.45 

1.00 

45.19 

0.0050 

0.0107 

0.627 


160 

3.30 

1.00 

103.07 

0.0311 

0.0454 

0.796 

0.529 

173 

2.87 

0.96 

159.73 

0.0326 

0.0365 

0.775 

0.674 

196 

3.41 

0.95 

188.30 

0.0345 

0.0381 

0.805 

0.703 

219 

3.89 

0.94 

170.32 

0.0385 

0.0459 

0.816 

0.659 

233 

3.83 

0.97 

148.91 

0.0417 

0.0472 

0.872 

0.601 

247 

2.73 

0.99 

121.08 

0.0450 

0.0548 

0.853 

0.497 

285 

0.10 

0.86 

2.55 

0.0005 

0.0004 

0.650 

0.034 
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Table 4 , Annual averages ( 152< DOY< 259) of the median grid cell values of LAI, CF, 
photosynthetic pigment concentration, FAPAR canopy ,FAPAR C hi, LUE ccmopy and LUE c ),i for the area 
around the Old Aspen flux site. 



LAI 

CF 

Photosynthetic 

pigment 

(pg/cm 2 ) 

FAP AR C anopy 
0-1 

FAPARci 

0-1 

LUE C anopy 
(ji mol C 
1 1 mol ’PPFD) 

LUE C hi 
(|i mol C 
|4 mol 1 PPFD) 


3.49 

0.98 

175.73 

0.856 

0.652 

0.0183 

0.0242 

2002 

3.35 

0.96 

161.32 

0.815 

0.617 

0.0183 

0.0245 

2003 

3.23 

0.94 

123.68 

0.794 

0.582 

0.0163 

0.0225 

2004 

3.54 

0.96 

140.32 

0.811 

0.623 

0.0244 

0.0310 

2005 

3.14 

0.97 

135.56 

0.820 

0.589 

0.0192 

0.0267 
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Figure Captions 

Fig. 1. Nadir MODIS NIR| reflectances (VZA < 5°) for the 5x5 area (2.5 km x 2.5 km) on 
DOY 224 in 2001. The central grid cell covers the Southern Old Aspen [SOA] tower site in 
Canada. 

Fig. 2. The relationship of APAR ch i = [(90 min PPFD)*FAPAR ch i] to APAR canopy = [(90 min 
PPFD) * F AP AR cai ,o Py | for 2001-2005. Solid diamonds (♦) indicate values retrieved during the 
core growing season (DOY:152-259; r = 0.87; APAR c hi = 0.713 *APAR can opy)- Early 
(DOY<152) and late (DOY>259) season values were combined and are indicated with open 
triangles (A) (r = 0.55; APAR ch i = 0.339* APAR can0 p y ). For all values (dashed line), r = 0.72 
and APAR C hi = 0.649*APAR canopy . 

Fig. 3. The relationship of GEP to APAR canopy = [(90 min PPFD)*FAPAR can0 py] in 2001-2005. 
Solid diamonds (♦) are data from during the core growing season (DOY 152-259). Open 
triangles (A) indicate values obtained before (DOY<152) or after (DOY>259). The apparent 
LUEcanopyi the slope of the relationship (0.0155 pmol C pmol' 1 PPFD, r 2 = 0.33), is similar to 
the core growing season value (0.0173) but is considerably higher than for days having low 
GEP in the spring and fall. 

Fig. 4. The relationship of GEP to APAR ch i = [(90 min PPFD)*FAPAR c hi] in 2001-2005. Solid 

diamonds (♦) are data from during the core growing season (DOY 152-259, GEP = 

0.0243*APAR ch i, r 2 = 0.63). Open triangles (A) indicate values obtained before (DOY<152) or 
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after (DOY>259) (GEP = 0.0208*APAR chl , r 2 = 0.64). The apparent LUE chl , i.e., the slope of the 
relationship for all values, is 0.0241 jimol C junol' 1 PPFD (r 2 = 0.78). 

Fig. 5. Comparison of the annual means ± SE for APAR C hi (♦) and APARcanopy (□) during the 
five-year period (2001-2005) for the growing season between DOY=l 52-259. APAR C hi was 
significantly lower than APAR can0 py in every year, averaging 235 jimol PPFD m'V 1 less than 
APAR can0 p y . APARtower is not included in this figure because unlike APAR c hi and APAR canop y, 
it was obtained under all sky conditions, i.e., both clear and cloudy. 

Fig. 6. Comparison of the annual means ± SE for MODlS-derived LUE c hi(>) and LUE can „ P y (□) 
during the five-year period (2001-2005) for the core growing season between DOY =152-259. 
LUEchi was significantly higher than LUE Cimopy in every year, averaging 0.007 jimol C jimol' 1 
PPFD higher. Annual LUE t0W er values (Krishnan et al., 2006) (A) agree well with our LUE c hi, 

falling within the SE range in 4 of 5 years. The maximum LUE (£ max ) for broadleaf forests 

used by the MODIS PSN model is showm as a horizontal dashed line. 

Fig. 7. Linear relationship between annual average LUE as a function of the cumulative 
precipitation for the core growing season (DOY= 152-259) in 2001-2005. 
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Fig. 1 . Nadir MODIS NIR 1 reflectances (VZA < 5°) for the 5 x 5 area (2.5 km x 2.5 km) on DOY 224 in 2001 . The 
central grid cell covers the Southern Old Aspen [SOA] tower site in Canada. 
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Fig. 2. The relationship of APAR chl = [(90 min PPFD)*FAPAR chl ] to APAR canopy = [(90 min PPFD)*FAPAR canopy ] for 
2001-2005. Solid diamonds (♦) indicate values retrieved during the core growing season (DOY:1 52-259; r = 0.87; 
APAR chl = 0.713 *APAR canopy ). Early (DOY<152) and late (DOY>259) season values were combined and are 
indicated with open triangles (A) (r = 0.55; APAR chl = 0.339* APAR canopy ). For all values (dashed line), r = 0.72 

and APAR nhl = 0.649*APAR_ nm/ . 
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Fig. 3. The relationship of GEP to APAR canopy = [(90 min PPFD)*FAPAR canopy ] in 2001-2005. Solid diamonds (♦) 
are data from during the core growing season (DOY 152-259,GEP=0.0173*APAR canopy ,r 2 =0.16). Open triangles 
(A) indicate values obtained before (DOY<152) or after (DOY>259)(GEP=0.0063*APAR canopy ,r 2 =0.18). The 
apparent LUE canopy , the slope of the relationship (0.0155 pmol C mmol-1 PPFD, r 2 = 0.33), is similar to the core 
growing season value (0.0173) but is considerably higher than for days having low GEP in the spring and fall. 
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Fig. 4. The relationship of GEP to APAR chl = [(90 min PPFD)*FAPAR chl ] in 2001-2005. Solid diamonds (♦) are 
data from during the core growing season (DOY 152-259, GEP = 0.0243*APAR chl , r 2 = 0.63). Open triangles (A) 
indicate values obtained before (DOY<152) or after (DOY >259) (GEP = 0.0208*APAR chl , r 2 = 0.64). The apparent 
LUE chl , i.e., the slope of the relationship for all values, is 0.0241 pmol C pmol-1 PPFD (r 2 = 0.78). 
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Fig. 5. Comparison of the annual means ± SE for APAR chl (♦) and APAR canopy (□) during the five-year period 
(2001-2005) for the growing season between DOY=1 52-259. APAR chl was significantly lower than APAR canopy in 
every year, averaging 235 jimol PPFD rrr 2 s- 1 less than APAR canopy . APAR tower is not included in this figure 
because unlike APAR chl and APAR canopy , it was obtained under all sky conditions, i.e., both clear and cloudy. 
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Fig. 6. Comparison of the annual means ± SE for MODIS-derived LUE ch , (♦) and LUE canopy (□) during the five-year 
period (2001-2005) for the core growing season between DOY =152-259. LUE chl was significantly higher than 
LUE canopy in every year, averaging 0.007 pmol C pmoM PPFD higher. Annual LUE tower values (Krishnan etal., 

2006) (A) agree well with our LUE chl , falling within the SE range in 4 of 5 years. The maximum LUE (£ max ) used by 
the MODIS-PSN model is shown as a horizontal dashed line(— ). 





LUEchi or LUEcanopy ( pmol C (imor 1 PPFD) 


0.032 
0,030 
0.028 
0.026 
0.024 
0.022 
0.020 
0.018 
0.016 

100 150 200 250 300 350 400 450 

Precipitation (mm) 



Fig. 7 Linear relationship between annual average LUE as a function of the cumulative precipitation for the core 
growing season (DOY=1 52-259) in 2001-2005. 




